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Abstract 

As  part  of  the  currently  postponed  National 
Wind  Tunnel  Complex  (NWTC)  program,  a  compu¬ 
tational  fluid  dynamic  (CFD)  study  of  a  large  trans¬ 
port  configuration  was  performed  using  the  AEDC 
chimera-overset  grid  XAIR  Euler  code.  Wall  poros¬ 
ity,  porosity  control  for  the  reduction  of  wall  interfer¬ 
ence,  correction  approach,  and  the  impact  of  refer¬ 
ence  point  selection  on  wall  interference  at  tran¬ 
sonic  conditions  (Mach  Number  0.85)  were  studied. 
These  effects  were  computationaliy  investigated 
using  a  clean-wing  model  representative  of  the  MD- 
1 1  aircraft  which  spanned  80  percent  of  the  test 
section  width.  Global  wall  corrections  were  deter¬ 
mined  by  a  streamtube  calculation  of  the  interfer¬ 
ence  flow  field.  Results  of  these  computations 
showed  that  global  corrections  for  Mach  number  at 
the  1/4  chord,  50-percent  semi-span  location  with 
angie-of-attack  correction  based  on  matching  wing 
lift  gave  very  low  interference.  Investigation  of  con¬ 
trol  of  porosity  showed  that  very  few  zones  were 
required  to  achieve  superior  results  for  global  cor¬ 
rections  and  for  model  pressure  distribution.  Resid- 
uai  blockage  gradient  and  flow-curvature  effects 
are  shown  to  be  treatable  by  very  small  changes  in 
side-wall  divergence  angle  and  model  geometric 
changes  that  are  equivalent  to  shearing  the  coordi¬ 
nate  system  aft  of  the  wing-fuseiage  juncture. 

Nomenclature 

Ca  Axial-force  coefficient 
Cd  Drag  coefficient 
Cl  Lift  coefficient 


Pitching-moment  coefficient 
Cn  Normal-force  coefficient 

Cp  Pressure  coefficient 

M  Mach  number 

R  Wali  resistance  [Eq.  (2)] 
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Subscripts 

T  tunnel  condition 

«>  free-air  condition 

1.0  Introduction 

The  facilities  for  the  National  Wind  Tunnel  Com¬ 
plex  (NWTC)  were  being  designed  with  the  expec¬ 
tation  that  the  maximum  design  Reynolds  number 
would  be  achieved  by  testing  models  with  spans 
up  to  80  percent  of  the  applicable  test-section 
dimension.  A  consequence  of  the  relatively  large 
model-to-tunnel  size  ratio  is  that  interference  pro- 
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duced  by  the  wind  tunnel  walls  may  be  unaccept¬ 
ably  large.  The  80  percent  span  ratio  is  believed  to 
be  the  largest  value  that  can  be  tested  with  any 
reasonable  chance  of  obtaining  high-quality  data. 
Data  quality  was  of  paramount  importance  for  the 
NWTC;  therefore,  wall  correction  and  wall  control 
strategies  were  being  investigated  as  an  integral 
part  of  the  development  process.  The  preliminary 
criteria  for  wall  interference  at  Mach  number  0.85 
were  established  for  the  uncorrected  levels  of 
blockage  (AM  <  0.003),  angle  of  attack  (Aa  ^  0.15 
deg  at  Cm  =  1 .0),  force  coefficients  (ACfg  <  0.005 
and  ACa  ^  0.0003),  and  pitching-moment  coeffi¬ 
cient  corrections  (AC,n  <  0.01).  The  intent  of  these 
criteria  was  to  have  the  level  of  wall  interference 
sufficiently  low  that  reliable  corrections  could  be 
made  to  the  data  to  support  meeting  the  overall 
uncertainty  in  the  measurement  of  one  drag  count, 
(Cd  =  0.0001).  The  purpose  of  the  CFD  evalua¬ 
tions  of  wall  interference  was  to  support  validation 
of  both  the  reasonableness  of  the  interference  and 
the  development  of  the  design  features  of  the  test 
section.  The  current  paper  summarizes  the  results 
obtained  and  documented  in  Ref.  1 . 


All  computations  were  solved  for  flows  about  a 
generic  representation  of  an  MD-11  clean-wing 
configuration.  All  in-tunnel  calculations  were  per¬ 
formed  for  a  model  that  is  scaled  to  span  80  per¬ 
cent  of  the  tunnel  width.  The  test-section  geometry 
of  this  tunnel  was  16  ft  wide  by  13  ft  high  with  no 
corner  fillets.  The  scaled  model  was  13.4  ft  in 
length,  and  the  model  nose  was  positioned  on  cen¬ 
terline  16.3  ft  downstream  of  the  test-section 
entrance.  The  test  section  was  42  ft  long.  The 
scaled  full-span  configuration  resulted  in  a  solid 
blockage  ratio  of  approximately  0.60  percent. 

For  all  flow  conditions  under  consideration,  lat¬ 
eral  symmetry  was  assumed  about  the  center  plane 
of  the  model.  Therefore,  only  half  of  the  flow  field 
was  modeled.  Figure  1  shows  a  portion  of  the  near¬ 
field  grid  system  for  the  MD-1 1 .  The  entire  system 
contains  nine  grids  with  1.1x10®  grid  points. 

3.0  Wall  Interference  Procedure 

In  the  current  investigation,  wall  interference 
effects  were  predicted  by  two  different  methods, 
constant  conditions  and  global  corrections. 


2.0  Computational  Procedures  3.1  Method  1 :  Constant  Conditions 


Wall  interference  estimates  were  based  on 
Euler  CFD  flow-field  solutions  which  were  per¬ 
formed  using  the  AEDC  chimera-overset  grid  XAIR 
code.^"^ 


In  the  constant  conditions  method,  the  wall 
interference  effects  are  predicted  at  a  constant 
geometric  angle  of  attack  and  upstream  Mach 
number  by  taking  the  difference  between  two  CFD 


Fig.  1 .  Chimera  inviscid  grid  system. 
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flow-field  solutions.  The  first  solution  is  a  free-air 
flow-field  calculation  where  the  outer  boundary  of 
the  computational  domain  extends  well  beyond  the 
location  of  the  tunnel  walls.  At  the  outer  boundary 
of  this  domain,  free-stream  conditions  are  pre¬ 
scribed.  The  upstream  and  downstream  outer 
boundaries  are  approximately  four  body  lengths 
upstream  and  two  body  lengths  downstream  of  the 
model,  respectively.  The  vertical  and  horizontal 
planes  of  the  outer  boundary  are  approximately 
five  tunnel  heights  from  the  centerline  axis  of  the 
model.  The  second  solution  is  an  in-tunnel 
flow-field  calculation  at  the  same  flow  conditions 
where  a  wall  boundary  condition  (Sec.  4.0)  is  pre¬ 
scribed  at  the  tunnel  walls.  Here,  the  flow  at  the 
test  section  entrance  is  prescribed  as  uniform  free- 
stream  flow,  while  the  flow  at  the  test  section  exit  is 
specified  by  an  extrapolated  boundary  condition.  A 
pair  of  flow-field  solutions  must  be  computed  and 
compared  to  determine  the  wail  interference  at 
each  selected  flow  condition.  To  minimize  any  pos¬ 
sible  difference  due  to  grid  resolution,  meshes 
internal  to  the  tunnel  are  identical.  The  local  effects 
of  wall  interference  on  the  model  are  examined  by 
subtracting  the  calcuiated  local  pressure  coefficient 
for  the  in-tunnel  solution(Cpx)  from  the  corre¬ 
sponding  free-air  solution(Cp  J  as 

ACp  =  Cp„-Cpj  (1) 

and  the  global  effects  on  forces  and  moments 
directly  follow  by  pressure  integration  of  the  individ¬ 
ual  solutions  and  subtracting  the  result. 

It  is  important  to  note  that  the  wali  interference 
procedure  outlined  above  determines  corrections 
to  the  model  forces  and  moments  at  the  fixed  tun- 
nei  free-stream  flow  conditions  and  at  a  fixed  geo¬ 
metric  angle  of  attack.  The  result  from  this  compu¬ 
tation  can  be  compared  to  the  NWTC  requirement 
for  acceptable  level  of  wall  interference.  Difference 
in  viscous  flow  development  obviously  will  affect 
the  results.  However,  the  method  could  include  an 
integral  boundary-layer  technique  to  make  a  vis¬ 
cous  correction 

3.2  Method  2:  Global  Corrections 

A  more  widely  used  wall  interference  procedure 
is  to  globally  correct  for  blockage  and  flow-curva¬ 


ture  effects  of  the  tunnel  walls  by  determining 
adjusted  remote  flow  conditions.  The  objective  of 
this  procedure  is  to  find  equivalent  free-aIr  flow  con¬ 
ditions  that  best  match  the  tunnel  data  or,  in  this 
case,  the  in-tunnel  calculations.  The  adjusted, 
equivalent  conditions  are  typically  determined  by 
selecting  a  blockage  correction  (AM)  and  flow- 
angle  correction  (Aa)  to  the  free-air  remote  condi¬ 
tions  which  force  the  free-air  and  in-tunnel  calcula¬ 
tions  to  match  at  a  given  reference  point  (or  points) 
on  the  model.  In  the  present  context,  this  procedure 
involves  a  minimum  of  two  additional  calculations. 

The  first  additional  flow-field  calculation  deter¬ 
mines  the  adjustments  to  angle  of  attack  and  Mach 
number  (AM  and  Aa).  This  is  done  by  performing 
an  inverse  computation  of  the  interference  flow 
field.  The  inverse  computation  involves  a  tunnel- 
empty  flow-field  or  streamtube  solution  with  the 
interference  pressure  distribution  specified  as  the 
boundary  condition  at  the  tunnel  boundaries.  This 
interference  pressure  distribution  is  determined  by 
subtracting  the  free-air  pressure  distribution  from 
the  tunnel  pressure  distribution,  both  of  which  are 
calculated  at  the  tunnel  walls  from  the  aforemen¬ 
tioned  free-air  and  in-tunnel  calculations.  The 
resulting  solution  of  the  inverse  flow  field  is  interpo¬ 
lated  onto  the  phantom  model  surface  to  provide 
disturbances  (AM  and  Aa)  at  selected  reference 
points.  A  complete  parametric  analysis  of  the 
selection  of  the  reference  point  location  is  beyond 
the  scope  of  this  investigation.  Typically,  the  flow 
corrections  are  determined  at  a  given  spanwise 
location.  At  the  given  spanwise  location,  a  different 
reference  point  is  selected  for  determining  AM  and 
Aa.  The  reference  points  are  selected  as  the  1/4- 
chord  location  for  AM  and  the  3/4-chord  location  for 
Aa.  The  selection  of  these  chordwise  locations  is 
based  on  classical  theory.^  The  3/4  chord  for  Aa  is 
the  classical  point  of  choice  to  minimize  pitching- 
moment  residuals.  Transonic  flow  features  may 
require  some  other  choice.  Two  spanwise  loca¬ 
tions,  25-percent  and  50-percent  semispan,  were 
selected  to  investigate  different  correction  strate¬ 
gies.  Initial  computations  were  performed  using  the 
corrections  that  were  determined  from  both  span- 
wise  locations.  As  wili  be  shown  in  the  results,  the 
50-percent  location  yielded  superior  blockage  cor¬ 
rections  and  was  selected  as  the  location  of  choice 
for  the  remaining  computations. 
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The  second  additional  flow-field  calculation  is 
performed  to  determine  the  free-air  flow  field  at  the 
adjusted  remote  conditions  (M  +  AM  and  a  +  Aa). 
Using  Eq.  (1),  the  residual  interference  is  deter¬ 
mined  by  comparing  the  free-air  flow-field  calcula¬ 
tion  at  adjusted  conditions  and  the  in-tunnel  calcu¬ 
lation  at  the  original  conditions  (appropriately  scal¬ 
ing  load  and  pressure  coefficients  for  the  difference 
in  dynamic  pressure  and  adjusting  stability-axis 
forces  for  differences  in  angle  of  attack). 

The  process  to  determine  the  optimum  remote 
conditions  is  an  iterative  procedure  where  the 
inverse  and  free-air  calculations  are  repeated  until 
the  residual  interference  is  minimized.  In  the  cur¬ 
rent  investigations,  only  one  inverse  computation 
was  required  to  obtain  an  accurate  blockage  or 
Mach  number  correction.  Additional  adjustments 
were  made  to  the  remote  angle  of  attack  to  best 
match  the  wing  lift  between  the  in-tunnel  and  free- 
air  calculations.  This  is  similar  to  the  technique 
used  by  Rizk.®  The  lift  matching  process  typically 
took  one  additional  free-air  computation.  However, 
because  the  interference  on  the  model  is  spatially 
variant,  residual  wall  interference  still  exists  in  the 
lift  contribution  from  the  other  model  components 
and  in  the  other  force  and  moment  coefficients  at 
these  adjusted  flow  conditions.  The  residual  is  pre¬ 
sumed  to  be  significantly  less  than  the  total  correc¬ 
tion  determined  for  the  original,  unadjusted  flow 
conditions. 

One  advantage  of  this  procedure  is  that  differ¬ 
ences  in  pressure  gradient  typically  are  much  less 
than  for  Method  1.  Hence,  incremental  corrections 
to  account  for  viscous  flow  development  differ¬ 
ences  are  expected  to  be  substantially  smaller. 

4.0  Wall  Boundary  Conditions 

The  porous  slotted  wall  behavior  was  simulated 
by  assuming  and  imposing  a  linear-homogeneous 
boundary  condition  representation  at  the  test  sec¬ 
tion  walls.  Provided  that  a  slotted  wall  has  suffi¬ 
cient  slots,  the  local  effect  of  the  discrete  slots  rap¬ 
idly  diminishes,  and  the  behavior  of  the  wall  is 
essentially  uniform  and  continuous  at  a  short  dis¬ 
tance  from  walls.^  The  NWTC  slotted  wall  configu¬ 
ration  was  designed  to  have  sufficient  slots  and  to 
behave  in  this  manner.  Therefore,  homogeneity 


was  a  reasonable  assumption  for  this  study.  A 
boundary  condition  was  selected  that  has  been 
shown  to  be  representative  of  the  wall  behavior  in 
the  NASA/ARC  1 1-fl  Tunnel  and  was  used  for  in¬ 
tunnel  calculations.  For  this  boundary  condition, 
the  cross-flow  behavior  at  the  wall  is  presaibed  as 


where  R  is  the  wall  resistance  or  porosity  parame¬ 
ter,  T  is  wall  porosity  or  wall  open  area  ratio,  and  0 
is  the  local  flow  angle  at  the  wall.  For  all  in-tunnel 
flow-field  calculations,  R  =  19.  To  simulate  test  sec¬ 
tion  walls  with  global  uniform  porosity,  the  value  of 
T  was  constant  over  all  walls.  To  simulate  variable 
porosity,  x  was  varied  locally  on  all  walls.  This 
boundary  condition  was  incorporated  into  the 
AEDC  XAIR/chimera  code  and  implemented  by  two 
methods.  The  first  method  calculates  0  at  each  iter¬ 
ation  in  the  flow-field  solution,  determines  the  wall 
pressure  from  Eq.  (2),  and  then  updates  and 
imposes  the  internal  energy  as  the  boundary  condi¬ 
tion  (along  with  extrapolated  values  of  density  and 
momentum).  The  second  method  calculates  Cp, 
determines  0  from  the  flow-field  solution,  and 
updates  and  specifies  the  appropriate  component 
of  momentum  as  the  boundary  condition.  The  sec¬ 
ond  method  was  utilized  when  variable  porosity 
was  investigated,  while  the  first  method  was  used 
for  all  uniform  porosity  computations. 

5.0  Results 

The  in-tunnel  and  free-air  flow-field  calculations 
about  the  model  were  performed  on  the  AEDC 
Convex  3800.  These  inviscid  computations  typi¬ 
cally  required  12  CPU  hours  to  converge.  Steady- 
state  convergence  was  considered  met  when  all 
calculated  model  forces  and  moments  coefficients 
remained  unchanged  to  five  significant  digits. 

5.1  Wall-Interference  Computations  at 
Constant  Conditions  (Method  1) 

Baseline  wall  interference  computations  for 
several  wall  configurations  and  flow  conditions 
were  performed.  This  task  involved  making  a  num¬ 
ber  of  in-tunnel  flow-field  calculations  at  different 
globally  uniform,  fixed  values  of  wall  porosity  for 
transonic  speed  conditions  and  the  corresponding 
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free-air  flow-field  calculations.  The  purpose  of 
these  calculations  was  to  quantify  the  level  of  wall 
interference  that  can  be  expected  for  walls  with 
globally  uniform  porosity.  All  calculations  were  per¬ 
formed,  as  well  as  the  determination  of  wall  inter¬ 
ference  increments,  at  the  fixed  transonic-speed 
conditions  of  =  0.85,  a  =  0  and  4  deg.  In  addi¬ 
tion,  in-tunnel  computations  were  made  at  discrete 
porosities  of  t  =  4 , 6  ,  and  9  percent. 


Table  1  shows  the  baseline  computations,  the 
flow  conditions,  and  the  resulting  forces  and 
moments  wall  interference  increments  using 
Method  1 .  The  comparisons  between  the  in-tunnel 
calculations  and  the  corresponding  free-air  model 

Table  1 .  Wall  Interference  Increments  for 

Baseline  Calculations  (Constant 
Conditions  Method) 
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Fig.  2.  Model  pressures  at  fixed  remote  condi¬ 
tions,  M  =  0.85,  a  =  4  deg,  t  =  4  percent, 
globally  uniform. 


pressure  distributions  for  a  =  4  deg  are  shown  in 
Figs.  2-4.  The  local  pressure  distributions  show  a 
significant  longitudinal  and  spanwise  variation  in 
the  increment  between  the  free-air  and  in-tunnel 
calculations.  The  wing-pressure  differences 
between  the  in-tunnel  and  free-air  calculations 
appear  to  be  minimized  at  4  percent  (Fig.  2).  As  the 
walls  are  opened  for  the  transonic  conditions,  the 
effect  on  the  model  wing  pressure  distribution  can 
be  seen.  The  pressure  level  increases,  the  shock 
strength  decreases,  and  the  shock  moves  fonvard. 

Figure  5  summarizes  the  calculated  wall  inter¬ 
ference  force  and  moment  increments  for  all  the 
inviscid  baseline  calculations  and  shows  the  varia¬ 
tion  of  those  increments  with  porosity.  The  base¬ 
line  calculations  show  that  there  is  a  significant 
variation  in  force  and  pitching-moment  increments 
with  uniform  porosity  and  that  wail  interference  var¬ 
ies  significantly  with  model  attitude.  Because  of  the 
gradient  of  wall  interference  over  the  model,  mini¬ 
mizing  all  forces  and  moments  simultaneously  can¬ 
not  be  accomplished  with  a  uniform  porosity.  With 
uniform  porosity  and  no  corrections  for  blockage  or 


Fig.  3.  Model  pressures  at  fixed  remote  condi¬ 
tions,  M  =  0.85,  a  =  4  deg,  x  =  6  percent, 
globally  uniform. 
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Fig.  4.  Model  pressures  at  fixed  remote  condi¬ 
tions,  M  =  0.85,  a  =  4  deg,  x  =  9  percent, 
globally  uniform. 


flow  curvature,  the  interference  increments  do  not 
meet  the  NWTC  wall  interference  criteria  for  tran¬ 
sonic  flow  conditions. 

Controllable  walls  were  considered  as  a 
method  to  meet  the  stringent  NWTC  data  quality 
specifications.  As  seen  in  the  baseline  computa¬ 
tions,  uniform  porosity  does  not  provide  the  desired 
data  quality  for  transonic  flows.  Variable  porosity 
was  investigated  to  show  the  possible  benefits  of 
controllable  walls  in  meeting  these  data  quality 
specifications.  Two  variable-porosity  methods 
were  computationally  investigated  using  the  linear- 
homogeneous  boundary  conditions.  In  the  first 
method,  it  was  assumed  that  there  is  an  infinite 
number  of  control  segments  over  the  walls.  With 
this  method,  the  porosity  varies  locally  over  the  test 
section  walls  to  best  match  the  local  wall  behavior 
determined  from  the  free-air  computation.  This 
method  will  be  referred  to  as  distributed  porosity 
control  and  represents  the  absolute  best  that  can 
be  done  with  variable  porosity  alone.  In  the  second 
method,  limited  control  was  imposed.  The  wall  was 
segmented,  and  the  porosity  was  fixed  laterally  or 
vertically  over  the  patch  and  allowed  to  vary  lin¬ 
early  in  the  longitudinal  direction.  The  longitudinal 
variation  was  constrained  to  be  piecewise  continu- 


m  (C  m  oo — C  fnT ) 


Fig.  5.  Wall  interference  increments  at  fixed  remote  conditions. 
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ous  between  patches.  This  method  will  be  referred 
to  as  segmented  porosity  control  and  represents  a 
more  realistic  implementation  of  variable  wall  con¬ 
trol.  Both  of  these  were  utilized  to  investigate  vari¬ 
able  porosity  for  the  two  transonic  flow  cases. 

From  the  free-air  baseline  calculations,  ideal 
values  of  the  local  porosity  for  the  transonic  cases 
were  determined.  This  was  accomplished  by  inter¬ 
polating  the  free-air  flow  field  onto  a  surface  equiv¬ 
alent  to  the  tunnel  boundary  and  then  calculating 
the  local  pressure  and  flow  angle  distributions  to 
determine  the  ideal  local  cross-flow  properties. 
Using  Eq.  (2),  the  local  values  of  porosity  were 
determined.  Practical  considerations  (e.g.,  the 
porosity  cannot  be  negative)  made  it  necessary  to 
constrain  the  porosity  such  that  1 0%  S  t  >  0%  (solid 
wall).  This  distribution  of  porosity  was  imposed  as 
the  boundary  conditions  and  the  in-tunnel  flow-field 
solutions  were  calculated.  The  results  of  these  cal¬ 
culations  were  compared  to  the  corresponding  free- 
air  solutions.  The  force  increments  are  shown  in 
Table  2.  The  wall  interference  is  significantly 
reduced  over  the  uniform  porosity  baseline  cases, 
as  can  be  seen  in  Fig.  5.  The  model  pressure  com¬ 
parisons  for  a  =  4  deg  are  shown  in  Fig.  6.  The 
results  from  these  calculations  show  a  significant 
reduction  in  wall  interference  by  incorporating  vari¬ 
able  porosity  control.  Figure  7  shows  the  corre¬ 
sponding  wall  cross-flow  distribution  (flow  angle  vs. 
pressure).  The  wall  cross-flow  distribution  for  the 
distributed  porosity  cases  and  the  free  air  shows 
good  agreement  except  for  a  slight  overshoot  in 
pressure  coefficient  at  the  model  center.  This  over¬ 
shoot  is  caused  by  an  acceleration  of  the  flow  from 
the  upstream  solid  wall  patch.  The  results  of  the  dis¬ 
tributed  porosity  computations  show  a  significant 
benefit  to  utilizing  variable  porosity. 

A  segmented  controllable  wall  was  investigated 
by  developing  a  patchwork  of  controllable  seg- 


Tabie  2.  Wall  Interference  Increments  for  Variable 
Porosity  Calculations  (Constant  Correc¬ 
tions  Method) 


Moo 

Mr 

If 

ay, 

deg 

percent 

ACp 

ACl 

0.85 

0.85 

0 

0 

dist 

-0.0008 

-0.0010 

0.0038 

0.85 

0.85 

0 

0 

seg 

-0.0008 

0,0018 

0.0023 

0.85 

0.85 

4 

4 

dist 

-0.0010 

0.0019 

0.0067 

0.85 

0.85 

4 

4 

seg 

-0.0008 

0.0040 

0.0064 

dist  =  distributed  porosity  variation 
seg  =  segmented  porosity  variation 


X/CorX/L 

Fig.  6.  Model  pressures  at  fixed  remote  condi¬ 
tions,  M  =  0.85,  a  =  4  deg,  distributed 
porosity. 
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ments.  The  goal  was  to  simulate  a  wall  where  the 
flow  through  the  slots  is  throttled  by  using  movable, 
plenum-side  cut-off  plates.  The  cut-off  plates  would 
allow  the  porosity  to  be  linearly  varied  in  the  flow 
direction.  In  addition,  the  goal  was  to  develop  a  seg¬ 
mented  wall  that  could  be  feasibly  implemented. 
The  patchwork  was  developed  by  examining  the 
wall  distributions  and  decid¬ 
ing  where  the  most  and 
least  wall  control  is  needed. 

The  most  control  is  needed 
in  regions  of  highest  gradi¬ 
ents,  e.g.,  model  wing, 
while  the  least  control  is  in 
the  regions  of  minimum  gra¬ 
dients,  e.g.,  upstream  and 
downstream.  Figure  8 
shows  the  patchwork  used 
for  the  computations.  This 
configuration  shows  six 
streamwise-controllable 
sections  on  the  top  and  side 
walls  and  one  on  the  bottom 
wall.  Implementation  of  this 
configuration  would 
require  six  controllable 
segments  along  all  or  selected  slots  on  the  top  and 
side  walls  and  one  control  segment  for  all  or 
selected  bottom  wall  slots. 

The  longitudinal  porosity  variation  was  deter¬ 
mined  by  anchoring  the  porosity  at  the  upstream 
and  downstream  ends  to  be  equal  to  the  interpo¬ 
lated  free-air  porosity  and  linearly  interpolating  over 
the  length  of  the  segment.  This  constrains  the  dis¬ 
tribution  to  be  piecewise  continuous  between  seg¬ 
ments.  The  wall  interference  results  of  these  calcu¬ 
lations  are  also  shown  in  Table  2  as  well  as  in  Fig 
5.  The  wall  interference  increments  are  comparable 
to  the  ideal  distributed  porosity  results.  The  model 
pressure  comparison  in  Fig.  9  and  corresponding 
wall  signature  in  Fig.  10  are  almost  identical  to  the 
distributed  porosity  results.  These  results  show  that 
limited  segmented  control  can  reduce  wall  interfer¬ 
ence  as  well  as  infinite  control,  provided  the  seg¬ 
mentation  is  strategically  placed.  Additional  work 
would  be  required  to  determine  the  optimum  place¬ 
ment  and  minimum  number  of  control  segments, 
the  solution  to  the  pressure  overshoot  in  the  wall 
signature,  and  an  active  wall  control  strategy. 


5,2  Wall  Interference  Computations  with  Global 
Corrections  (Method  2) 

For  selected  uniform  and  segmented  porosity 
calculations,  adjustments  to  free-air  flow-field  angle 
of  attack  and  Mach  number  were  determined  to  find 
equivalent  free-air  conditions  that  best  match  the 


Fig.  9.  Model  pressures  at  fixed  remote  conditions, 
M  =  0.85,  a  =  4  deg,  segmented  porosity. 
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in-tunnel  calculations.  The  purpose  of  these  com¬ 
putations  was  to  determine  the  residual  interfer¬ 
ence  after  the  computations  have  been  corrected 
for  first-order  blockage  and  flow-curvature  effects. 
The  flow-field  adjustments  were  determined  using 


S’  0.50 

1 - 1 - 1  1 - 1 - 1 - 1 - 1 - r*  ■■"-n 

Free-Air  - 

«  0.25 

t  =  6% - 

B  0 

Segmented  T 

^  -0.25 

S  -0.50 

lowier  wall 

^>.06  -0.04  -0.02  0  0.02  0.04 

Cn 


the  procedure  described  in  Sec.  3.0.  Five  cases 
were  selected  and  are  shown  in  Table  3.  These 
cases  include  two  porosity  distributions  at  a  =  0  deg 
(a  uniform  porosity,  t  =  6  percent,  and  the  seg¬ 
mented  porosity  distribution)  and  three  porosity  dis¬ 
tributions  at  a  =  4  deg  (two  uniform  porosities,  t  =  4 
and  6  percent,  and  the  segmented  distribution 
porosity). 

Figure  11  shows  representative  streamtube 
results  for  the  M  =  0.85,  a  =  4  deg,  and  x  =  6  per- 


Cp 


Fig.  10.  Wall  signature  comparison  at  M  = 
0.85,  a  =  4  deg,  segmented  porosity. 


cent  case  interpolated  onto  a  phantom  surface  of 
the  MD-1 1  wind/body/horizontal  tail  configuration. 
Shown  are  the  Mach  number  and  flow  curvature 
effects  of  the  walls.  The  contour  plots  show  a  sig¬ 
nificant  spanwise  and  longitudinal  gradient  in  the 

Table  3.  Wall  Interference  Increments  at  Adjusted  Free- 
Stream  Conditions  (Global  Corrections  Method) 


m 

a^,  deg 

Or, 

deg 

percent 

ACp 

ACl 

Ref. 

Loc., 

percent 

QJQgQ] 

■I 

6 

UlQIggl 

25 

■iina 

QESEl 

0 

6 

-0.0010 

-0.0023 

50 

1  0.85048 

BSMI 

0 

-0.0007 

-0.0006 

fiT;glrATil 

ESSI 

4 

6 

-0.0027 

-0.0103 

^■1 

EEaESI 

4 

6 

-0.0041 

0.0177 

50 

4 

4 

-0.0075 

50 

4 

seg 

l22sls!i£l 

50 

Tunnel  forces  and  moments  adjusted  to  the  dynamic  pressure  and  angle 
of  attack  of  the  free-air  conditions, 
seg  =  segmented  porosity  variation 


a.  Blockage  b.  Flow  angle  (radians) 

Fig.  1 1 .  Interference  flow-field  solution. 
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interference  field.  The  streamtube  results  were 
used  to  determine  the  adjustments  to  the  remote 
conditions. 

Initially,  remote  condition  corrections  were 
determined  from  two  spanwise  reference  locations 
in  the  streamtube  calculations.  These  spanwise 
locations  were  the  25-  and  50-percent  semispan 
locations.  Figures  12  and  13  show  the  model  pres¬ 
sure  agreement  for  these  initial  calculations,  while 
Table  3  shows  the  residual  force  and  moment  dif¬ 
ferences  at  the  corrected  conditions.  The  calcu¬ 
lated  in-tunnel  pressures  and  loads  were  all  cor¬ 
rected  to  the  free-air  dynamic  pressure  and  angle 
of  attack.  A  significant  improvement  in  the  wing 
pressure  comparison  can  be  seen,  particularly  for 
the  flow  determined  at  the  50-percent  semispan 
location  (compare  Fig.  13  to  Fig.  3).  The  pressure 
comparisons  clearly  show  the  advantage  of  using 
the  50-percent  semispan  location  as  the  reference 
point.  Similar  improvements  in  the  wing  pressure 
agreement  were  observed  for  the  other  fixed 
porosity  cases.  For  the  remaining  cases,  the  cor- 


Fig.  12.  Model  pressures  at  adjusted  remote  condi¬ 
tions,  M  =  0.847,  a  =  3.83  deg,  x  =  6 
percent,  globally  uniform,  reference  points 
located  at  25  percent  semispan. 


Fig.  1 3.  Model  pressures  at  adjusted  remote  condi¬ 
tions,  M  =  0.8485,  a  =  3.84  deg,  x  =  6 
percent,  globally  uniform,  reference  points 
located  at  50  percent  semispan. 

rections  to  the  remote  conditions  were  determined 
at  the  50-percent  semispan  reference. 

Model  pressure  comparisons  for  the  segmented 
porosity  are  shown  in  Fig.  14.  Pressure  agreement 
for  the  segmented  porosity  was  already  quite  good 
and  only  a  small  improvement  is  noticed. 

Although  a  significant  improvement  was  made 
in  the  wing  pressure  agreement,  pressure  differ¬ 
ences  still  exist  on  the  aft  fuselage.  The  differences 
are  primarily  due  to  the  longitudinal  gradients  in 
flow  curvature  and  velocity  between  the  wing  and 
the  tail.  These  gradients  can  be  seen  in  the  inverse 
computations  of  Fig.  1 1 . 

This  case  yielded  a  difference  in  flow  angle  of 
0.172  deg  between  the  wing  and  tail  reference 
points.  Normally,  the  difference  in  flow  angle  is 
resolved  by  modifying  the  tail  angle.  A  different 
approach  was  used  here.  The  whole  aft  fuselage 
was  warped  aft  of  the  wing  root.  Each  longitudinal 
grid  plane  downstream  of  the  wing  root  was  rotated 
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Fig.  14.  Model  pressures  at  adjusted  remote  condi¬ 
tions,  M  =  0.85096,  a  =  3.932  deg, 
segmented  x,  globally  uniform,  reference 
points  located  at  50  percent  semispan, 
about  a  wing  root  reference.  The  rotation  angle 
was  determined  by  a  linearly  weighted  flow  angle 
distribution  based  on  the  distance  from  the  wing 
root.  Figure  15  shows  the  model  pressure  compar¬ 
isons  for  the  warped  fuselage  calculations.  Signifi¬ 
cant  improvement  was  obtained  in  the  aft  fuselage 
agreement  (see  Fig.  13  for  comparison).  An 
improvement  was  also  seen  in  the  load  increments 
shown  in  Table  4.  The  process  of  warping  the  fuse¬ 
lage  is  similar  to  the  shearing  techniques  used  to 
account  for  aeroelastic  distortion. 

The  longitudinal  velocity  gradient  or  tunnel  buoy¬ 
ancy  is  normally  eliminated  by  adjusting  the 
wall  angle.  Side  wall  divergence  was  numer-  Table 
ically  attempted.  However,  problems  with 
coupling  the  side  wall  divergence  angle  and 


Force  and  moment  increments  at  the  given  con¬ 
ditions  are  compared  to  those  at  the  adjusted  con¬ 
ditions  in  Fig.  16.  The  load  increments  changed  sig¬ 
nificantly  for  the  adjusted  remote  conditions.  For  all 
cases  except  a  =  4  deg,  x  =  4  percent  and  a  =  0  deg, 
X  =  6  percent  (remained  unchanged),  the  magnitude 
of  the  lift  increment  decreased,  and  for  all  cases 
except  a  =  0  deg,  x  =  6  percent,  the  magnitude  of 
the**^pitching-moment  increments  increased.  The 
increasing  pitching  moment  increments  show  the 
need  to  include  tail  incidence  corrections  as  part  of 
the  correction  strategy.  With  a  combination  of  using 


Fig.  15.  Model  pressures  at  adjusted  remote  condi¬ 
tions,  M  =  0.8485,  a  =  3.84  deg,  x  =  6 
percent,  globally  uniform,  reference  points 
located  at  50  percent  semispan,  warped 
fuselage. 

Wall  Interference  Increments  for  Warped  Aft  Fuse¬ 
lage  and  Baseline  Configurations  (Global 
Corrections  Method) 


the  porous  wall  boundary  condition  were  not 
successfully  resolved  before  the  completion 
of  the  project.  The  correction  approach  to 
Mach  number  gradient  by  controlling  the  side 


wall  angle  is  practical  and  one-dimensional 
flow  analysis  revealed  only  a  small  change  to 
wall  angles  of  less  than  0.03  deg  is  needed. 


Tunnel  forces  and  moments  adjusted  to  the  dynamic  pressure  and  angle  of 
attack  of  the  free-air  conditions. 

BC  =  baseline  configuration 
WF  =  warped  fuselage 
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Fig.  16.  NWTC  wall  interference  comparison,  fixed  vs.  adjusted  conditions,  nominal 
M  =  0.285. 


Table  5.  Wall  Interference  Increments  for  ACL(wing)  Minimized 
(Global  Corrections  Method) 


0 

6 

25 

0.84900 

0 

6 

-0.0011 

-0.0018 

0.0081 

50 

0.85048 

-0.0326 

■■ 

seg 

-0.0007 

-0.0009 

0.0039 

50 

0.84700 

3.8960 

4 

6 

-0.0021 

-0.0037 

0.0195 

25 

0.84850 

3.8389 

4 

6 

-0.0014 

-0.0043 

0.0178 

50 

0.84940 

3.9452 

4 

4 

-0.0021 

-0.0053 

0.0229 

50 

4 

seg 

-0.0021 

ESStfil 

50 

Tunnel  forces  and  moments  adjusted  to  the  dynamic  pressure  and  angle  of 
attack  of  the  free-air  conditons. 
seg  -  segmented  porosity  variation 


variable  porosity,  side  wall  divergence  control,  and 
adequate  wall  corrections,  it  does  appear  that  the 
flow  quality  specifications  for  the  NWTC  can  be  met. 

Loads  were  calculated  for  each  model  compo¬ 
nent.  Forces  and  moments  contribution  for  the  wing, 
body,  horizontal  tail,  and  vertical  tail  were  tabulated 
for  each  run.  With  this  tabulation,  adjustments  to 
angle  of  attack  were  determined  to  minimize  the 
wing  lift  difference.  The  lift  difference  between  the 
in-tunnel  and  free-air  solution  was  driven  to  less 
than  one  count  in  lift  coefficient  (0.0001).  Once  the 
wing  lift  slope  was  determined  for  each  nominal 


angle  of  attack,  the  process  required  only 
one  additional  calculation.  The  purpose  of 
these  computations  was  to  explore  addi¬ 
tional  fine  tuning  of  the  above  correction 
strategy.  This  minimization  process  was 
performed  for  the  cases  in  Table  3  (as  well 
as  others).  Figures  17  and  18  show  the 
model  pressures  comparison.  The  resid¬ 
ual  interference  loads  increments  are 
shown  in  Table  5.  Comparison  of  the  force 
and  moment  increments  at  the  given  con¬ 
ditions  to  those  at  the  adjusted  conditions 
with  minimized  wing  lift  increments  is  shown  in  Fig. 
19.  No  significant  differences  are  seen  in  the  total 
force  and  moment  increments  between  these 
results  and  the  results  presented  in  Fig.  16.  Table  6 
shows  the  wing  force  and  moment  increments 
before  and  after  the  minimization  process.  Note  that 
the  wing  lift  is  matched. 

6.0  Conclusions 

A  computational  investigation  of  wall  interfer¬ 
ence  in  the  13  x  16-ft  NWTC  test  section  was  per¬ 
formed  for  transonic  conditions.  The  study  involved 
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Table  6.  Force  and  Pitching-Moment  Increments  Before  and  After  Wing-Lift  Minimization 
(Global  Corrections  Method) 


Conditions 

Before  Minimization 

Afl 

er  Mininu 

zation 

Ref. 

Mr 

percent 

deg 

aj, 

deg 

ACq 

ACl 

ttoc,  deg 

deg 

ACp 

AC„ 

ACl 

Loc., 

percent 

EES 

6 

tm 

B2SiSiiisi 

EE&SI 

ESE^I 

EEESS 

0 

25 

0.84900 

EES 

6 

-0.1000 

0 

-0.0003 

0,0017 

-0.0004 

o 

BESi 

0.0014 

0 

50 

0.85048 

EES 

seg 

— 

0 

-0.0001 

0,0006 

ESSEI 

-0.0326 

0 

QgggH 

0 

50 

0.84700 

6 

3.8300 

4 

-0.0014 

0.0089 

3.8960 

4 

qqqEI 

0 

25 

0.84850 

EES 

6 

3.8400 

4 

-0.0002 

0.0020 

0.0001  . 

3.8389 

4 

0.0001 

50 

0.84940 

EES 

4 

3.9182 

4 

0.0052 

3.9452 

EM 

-0.0007 

50 

0.85096 

EES 

seg 

3.9320 

4 

0.0000 

-0.0001 

0.0022 

3.9054 

0,0013 

0 

50 

Tunnel  forces  and  moments  adjusted  to  the  dynamic  pressure  and  angle  of  attack  of  the  free-air  conditions 
seg  =  segmented  porosity  variation 


Fig.  17.  Model  pressures  at  adjusted  remote  condi¬ 
tions,  M  =  0.8485,  a  =  3.8389  deg,  t  =  6 
percent,  globally  uniform,  ACl  (wing)  mini¬ 
mization  for  3.84  deg  case. 

computing  the  flow  over  a  clean-wing  MD-1 1  type 
model  that  spanned  80  percent  of  the  tunnel  width. 
Steady-state  converged  solutions  to  the  Euler 
equations  were  obtained  using  the  AEDC  XAIR 
code.  A  linear  homogeneous  wall  boundary  condi¬ 
tion  was  used  that  is  representative  of  the  slotted 
walls  in  the  NASA/ARC  1 1-ft  tunnel.  The  investiga¬ 
tion  involved  computing  a  number  of  baseline  cal¬ 
culations  to  determine  the  variation  of  wall  interfer¬ 
ence  with  porosity  at  fixed  geometric  angle  of 


Fig.  1 8.  Model  pressures  at  adjusted  remote  condi¬ 
tions,  M  =  0.85096,  a  =  3.9054  deg, 
segment  porosity,  globally  uniform,  ACl 
(wing)  minimization  for  3.932  deg  case. 

attack  and  Mach  number.  The  benefit  of  using  con¬ 
trollable  walls  was  investigated.  Wall  correction 
strategies  were  also  investigated. 

The  baseline  calculations  showed  that  there  is 
a  significant  variation  in  force  and  pitching-moment 
increments  with  uniform  porosity  and  that  wall 
interference  varies  significantly  with  model  attitude. 
Because  of  the  gradient  of  wall  interference  over 
the  model,  minimizing  all  forces  and  moments 
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Fig.  18.  NWTC  wall  interference  comparison,  fixed  vs.  adjusted,  minimized  ACl  (wing), 
nominal  M  =  0.85. 


simultaneously  cannot  be  accomplished  with  a  uni¬ 
form  porosity.  With  uniform  porosity  and  no  correc¬ 
tions  for  blockage  or  flow  curvature,  the  interfer¬ 
ence  increments  did  not  meet  the  wall  interference 
criteria  for  transonic  flow  conditions. 

Variable  porosity  schedules  based  on  the  free- 
air  solutions  were  determined  at  two  transonic  flow 
conditions.  Two  schedule  types  were  used,  one 
with  an  infinite  number  of  control  segments  and 
another  with  finite  control  segments.  The  results 
showed  that  variable  porosity  significantly  reduced 
the  wall  interference.  The  benefit  and  even  the 
necessity  of  using  controllable  walls  to  meet  data 
quality  specifications  was  clearly  shown.  Ideal  dis¬ 
tribution  of  porosity  with  infinite  number  of  control 
segments  yielded  comparable  results  with  limited 
finite  segmented  control.  Additional  work  is 
required  to  determine  the  optimum  placement  and 
minimum  number  of  control  segments,  the  solution 
to  the  pressure  overshoot  in  the  wall  signatures, 
and  an  on-line,  feed-fonward  wall  control  strategy 
to  predict  wall  settings.  The  current  investigation 
did  show  that,  with  a  limited  number  of  control  seg¬ 
ments,  wall  interference  can  be  greatly  reduced. 


Free-air  calculations  were  performed  at  adjusted 
remote  conditions  for  several  cases  to  examine  wall 
correction  strategies  and  to  determine  residual 
interference  after  first-order  Mach  number  and  flow 
curvature  effects  were  considered.  Adjusted  condi¬ 
tions  were  determined  at  given  reference  points 
from  an  inverse  computation  of  the  interference 
field.  Corrections  determined  at  two  different  span- 
wise  locations  showed  that  the  50-percent  semis¬ 
pan  location  yielded  superior  blockage  corrections. 
Wing  pressure  comparisons  at  the  adjusted  condi¬ 
tions  showed  a  significant  improvement  over  the 
pressure  comparisons  at  the  given  conditions.  For 
most  cases,  the  magnitude  of  the  lift  increment 
decreased,  but  the  magnitude  of  the  pitching- 
moment  increments  increased.  The  increasing 
pitching-moment  increments  showed  the  need  to 
include  tail  incidence  corrections  as  part  of  the  cor¬ 
rection  strategy.  The  tail  incidence  was  effectively 
changed  by  warping  the  aft  fuselage  to  account  for 
flow  curvature  gradients  between  the  wing  and  the 
tail.  Significant  improvement  was  obtained  in  the  aft 
fuselage  pressure  agreement  and  the  pitching- 
moment  increments.  Investigating  side  wall  control 
to  eliminate  tunnel  buoyancy  was  not  completed 
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and  should  be  performed  to  confirm  the  one-dimen¬ 
sional  analysis. 

With  a  combination  of  using  controllable  poros¬ 
ity,  possible  side-wall  angle,  and  adequate  wall 
corrections,  the  results  indicated  that  the  transonic 
wail  interference  specifications  for  the  NWTC  can 
be  met,  and  the  resulting  level  of  wall  corrections 
needed  is  small. 
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